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ABSTRACT 


An examination of a mathematical analysis of Pulse Tube refrig- 

 eration as proposed by Rea was carried out. Rea originally investi- 

gated the performance of thermal regenerators subjected to rapid 

pressure and flow cycling and found that this class of regenerators 

is characterized by two basic parameters; 

(1) The ratio of the convective heat transfer E the regenerator 
to the heat generated by compression, 

(2) The ratio of the rate of mass storage within the void volume of 
the regenerator to the through-flow. 

Rea's proposal was to treat the Pulse Tube as a very efficient regen- 

erator in series with a very inefficient regenerator, and, thus, 

characterize the performance of the Pulse Tube by the same two para- 

meters, 

The examination of Rea's proposal indicates that, in addition 
to the two parameters which describe the performance of an efficient 
regenerator, there is, for the case of the inefficient regenerator, a 
third parameter which depends upon the average pressure and the rate 
at which the pressure and mass flow cycle. 

An experimental evaluation of Rea's proposal modified to include 
the third parameter was performed on a Pulse Tube consisting of a 
spherical matrix regenerator in series with a void tube. There was 
good agreement observed between theoretical and experimental results. 


Thesis Supervisor: Joseph L, Smith, Jr. 


Title: Associate Professor of Mechanical Engineering 
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NOMENCLATURE 


cross-sectional area of the regenerator or pulse tube, 

ft", 

cross-sectional area of the regenerator or pulse tube open 
to flow (A, =€ A), #52. 

heat transfer surface area per foot of length of regenerator 
or pulse tube, ft^/ft. 

dimensionless variable defined by equation (A.l8). 
dimensionless constant defined by equation (2.31). 
integration constant m - Т. - му?” | 

constant defined by equation (2.19), lb, - орг, 

specific heat of regenerator matrix or pulse tube wall, 
BIU/1lb, - "F. 

constant pressure specific heat of the gas, BTU/1b, - Ор, 
constant volume specific heat of the gas, BTU/lb, = °F, 
gas enthalpy at a position x along the regenerator or pulse 
tube, BIU/lb.. 

heat transfer coefficient at position x along the regenera- 
tor or pulse tube, BTU/hr«ft? - F, 

heat transfer coefficient, BTU/hr-ft<-oF 。 

thermal conductivity, BIU/nr=ft-"F, 

regenerator or pulse tube length, ft. 

gas mass flow at position x along regenerator or pulse 


tube, lb,/hr. 


dimensionless mass flow (m š п/п). 


mix. 





m 


H 


HS 


r3 ct eU De 


mass flow at x = 0 in regenerator or pulse tube, 1b /hr. 
matrix mass per foot of regenerator length or mass per 
foot of pulse tube length, lb,/ft; | 
positive exponent whose value depends on regenerator 
matrix (n = 0 for laminar flow pulse tube). 

gas pressure in regenerator or pulse tube, 1br/ft<。 

heat flux, BTU/hr. 

gas constant, ft-1b,/1b,-"F. 

time, hours. 

temperature of the regenerator matrix or pulse tube wall 
at position x, oM 


temperature of the gas at position x along the regenera- 


tor or pulse tube, “F, 


internal energy of the gas, BTU/lb,. 

volume, ft’, 

dead volume at end of pulse tube, ft3 。 

volume of pulse tube, ft. 

void volume of regenerator, ft, 

total volume of the regenerator or pulse tube, ft^. 
dimensionless distence along regenerator or pulse tube 
defined by equation (2.26). 


distance along regenerator or pulse tube, ft. 


Greek Letters 


dimensionless Pulse Tube parameter defined by equation (3.18). 

heat transfer area per unit total volume of the regenerator 
AmL 

or pulse tube (P = =), ftc /ft3 。 


== К 





x 


D م‎ 


ANA 


-С 


Subscripts 


avg 


cond 


eff 
exper 
Б 
heater 
H 


losses 


dimensionless regenerator or pulse tube void fraction. 
dimensionless ratio of specific heats (Jv = с „су ). 
dimensionless pulse tube parameter defined by equation 
(2.35). 

gas density, lb/ft. 

dimensionless matrix or pulse tube wall temperature 
CO в). 

duration of cycle, hours. 

duration of a compression period, hours, 

duration of an expansion period, hours. 

percentage of the total cycle time devoted to compression 


or expansion, 


average, 

compression, 

conditions at cold end of Pulse Tube. 

conduction. 

expansion, 

effective, 

experimental, 

gas, 

used with Q to indicate heater input to Pulse Tube, 
conditions at hot end of pulse tube. 

used with Q to indicate conduction and radiation losses 
in Pulse Tube. 


refers to regenerator matrix or pulse tube wall. 


mat” 





max 


pt 


rad 


theor 


refers to maximum value of the variable w or p. 
refers to minimum value of the variable p. 

pulse tube or regenerator conditions at x - w * O. 
pulse tube, 

radiation, 

regenerator; also refers to refrigeration Cone 
theoretical. 


conditions at warm end of regenerator. 


Sa 





CARPTER І 
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A recent paper by Gifford and Longsworth (1)^ describes an inter- 
esting thermodynamic effect which frequently is ignored in prectice, 


E McC. ies a potentially useful application, As Gifford and 


FIGURE 1 


SCHzVATIC OF THE PULSE TUBE 


ال ^ ?^ 











Longsworth indicate (2), it is not generally recognized that the cycling 
of a gas in a closed volume can establish large temperature gradients 
within such a volume. Gifford and Longsworth, by suitably arranging a 


chermal regenerator, heat exchangers, and à void volume, have created a 


< 





« Numbers in ( ) refer to Refererces. 


M 





бөуісе which preserves these temperature gradients in an essentially 


3 ` 


Ah. е а А. 4242 a : x RAS „аа ` کہ‎ - 
вол с Бласе, oven though the sas is undergoing rapid cycling. 


U 


Pulse Tube. It is illustrated schematically 
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Tre construction and operation of the device are simple. “This 
simplicity when coupled with the vromising experimental results obtained 
tuus far, indicates that the Pulse Tube should find practical applica- 
ticn in tke near future, perhaps as a small cooler in electronic equiv- 
mente Gifford and Longsworth (2) report refrigeration across temperature 
intervals as larze as 2099 F. with a single-stage system. The units may 
also be oporated in multi-stages, so larger temperature differences ere 
entirely feasible. 

The Pulse Tube consists of a regenerator connected in series with 
a void tuse, which is the pulse tube,” The void volume of the regen- 


erator is Vas The volume of ine pulse tube is Vot” Gas enters and 


exits the cevice through the warm end oF the regenerator, which is at 
temperature Iye The volume Уш, which is connected to the end of the 
ша шсбе. is maintained at a constant temperature, Tue Tnis may be 
&ccompiished by circulating cooling water around the volume V.. 

The device operates in the following manner: tne inlet valve 


Spens acmitting high pressure gas into the regenerator, thereby 


inducing a mass flow into the regenerator and pulse tube. 





we Pulse Tube capitalized refers to the entire device (regenerator 
plus pulse tube); pulse tube without capital letters refers to 
tre void tube. Wherever the void tube is intended, but gram 
matically capital letters are required, the follcwing novation 
will be used: "pulse tube", 





concurrenuly), there is a continuous pressure rise throughout the Pulse 


& 


Juoe。 2% sore predetermined maximum pressure, ра within the device, 


E^ 


tne inet valve closes and the exit valve opens. This allows the gas 


within the Pulse Tube to flow back out through the warm end of ine regen- 
Колот. 

When the inlet valve is open, the gas in the pulse tube is heated 
as it is compressed by the entering gas. Some of this heated gas flows 


into the volume Vu? where it is cooled by tne cooling medium to ihe tem- 


nerature 1, When the exit valve opens, the gas is cooled further by 


1% 


«+ 


expansion as it flows back out of ihe volume V Some of this cool gas 


H° 
must precool the entering warm gas during the next compression half- 
cycle. The remainder of the cool gas is available for refrigeration. 


tion becomes available, once a steady cyclic state is reached, 


fo 


Reiriger 


at the connecting tubing between the regenerator and the pulse tube, 


-- 


9 
“nis refrigeration load is labeled ©з in Figure Y and is available at 


the temperature Tas 


Previous Work 





The first analysis of the pulse tube was given by Gifford and 
Longsworth in their first and second papers (1,2). This analysis did 
rot predict compatible pulse tube capacity results with those they had 
obtained experimentally. Gifford and Longsworth recognized that their 
original analysis of the pulse tube was inadequate. Tneir original 
analysis neglected any heat exchange between the gas and the pulse tube 
wall, They discovered that this heat exchange could contribute a sub- 
stantial amount to the refrigeration effect. In other words the pulse 


tube itself acts as an inefficient regenerator, The third paper by 





Gifford and Longsworth (3) discussed the aspects of this heat exchenge 
between the gas and the pulse tube wall in detail. However, no complete 
analysis including this effect was presented. The problem was primaril 


one of cetermining the appropriate temperature difference to use in 


FIGURE TI 
SCHEMATIC OF REA'S 2-PART MODEL FOR A REGENERATOR (L) 
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order to calculate correctly the heat transferred in the process under 
consideration. 

Concurrently with the work of Gifford and Longsworth, Rea was 
studying the performance of thermal regenerators subjected to rapid 
pressure and flow cycling (4). Reats objective was to develop a theory 


for predicting the axial temperature distribution, mass flow rate 


а 





dästrizution, and cffectiveness of thermal regenerators operated under 
the above conditions. Іп developing his theory Rea established his soe 
called 2=-Part Model (L), which is patterned after the lr ноп 
cycle (5), A schematic of Reats 2=Part Model is shown in Figure ITI. 

By making energy and mass balances for the gas and the matrix in 
crn elemental control volume, Rea was able to write the governing differ- 
ential equations for his problem, After applying appropriate simplifi- 
cations, Rea reduced his problem to a system of two differential equa- 
tions, which he ва numerically with the aid of a computer. A complete 
description of Rea!s 2-Part Model and the development of the governing 
differential equations are presented in Appendix A. 

When the work of Gifford and Longsworth became known to him, Rea 
recognized that the Pulse Tube operates in the same manner as that 
described by his 2=Part Model, Rea concluded that since a pulse tube 
was nothing more than an inefficient regenerator, his theory should 
apply equally well to the entire Pulse Tube as to a simple regenerator. 
Thus, Rea extended his theory to include the pulse tube. 

Rea tested his theory for the regenerator and got fairly accurate 
results. Rea, however, did not test his theory as extended to include 
the pulse tube. Thus, it is the purpose of this thesis to determine 
whether or not Rea's Theory correctly describes the performance of the 
Pulse Tube by, first, examining his theory and analysis for the entire 
Pulse Tube and by, second, measuring the refrigeration effect of an 


operating Pulse Tube. 





CHAPTER II 


EXAMINATION OF REA'S THEORY AS APPLIED 


TO A REGENFRATOR AND "pulse tube" 


Reats 2-Part Model (h) 

When faced with the task of formulating a mathematical analysis 
of a given system, the most important step is the selection of an 
appropriate model. Ideally one wants a model that is amenable to a 
theoretical treatment. On the other hand one does not want his model 
to be over-simplified, for then the problem may lose its significance. 
Rea selected as his model for the regenerator, and, hence, the pulse 
tube, a model that satisfies both of these criteria, This "2-Part 
Model", as Rea refers to it, is shown schematically in Figure II of 
Chapter I. The usefulness of this model is proved by the fairly 
accurate results which are obtained by its use and by the fact that 
certain dimensionless parameters which are definite aids in defining 
the problem evolve from its use, 

The operation of the 2-Pert Model is as follows: the inlet 
valve opens allowing high. pressure gas to enter at the top of the 
devices During this half-cycle the pressure within the device con- 
tinually increases with the mass flow rate out of the bottom of the 
device equal to P After a predetermined maximum pressure, P max, 
is reached within the device, the inlet valve closes and the exit 
valve opens, During the exhaust half-cycle the process described 
above is reversed, with mass flow rate into the bottom of the 


device equal to m. The pressure within the device decreases to 


“бө 





its minimum value, p during this half-cycle. After a period of 


min? 
time, a steady, cyclic condition is reached. 
In examining Rea's theory for the Pulse Tube, the theory for the 


individual components, ieee, the regenerator and pulse tube, will be 


examined separately, 


Везас Differential Equations for the Regenerator and "pulse tube! 


The governing differential equations are derived in detail in 
Appendix À by making a mass balance and energy balance for the gas end 
for the matrix in an elemental control volume, dx, at a position x 
within the regenerator and pulse tube respectively, The derivation of 
the regenerator equations is entirely attributable to Rea (lh). The 
modification to Rea's Theory as applied to a pulse tube is original. 


The basic equations for the regenerator and pulse tube are: 


Gas Energy Palance: 


= = (mT,) + ВАТ - Т) + со дР - 





р ed Cel 
9 mE 5 (2.1) 
Matrix Energy Balance: 
hA (T – Т) 2 Ме 91 
rT, ) Сп x+ 2 (2.2) 
Gas Continuity: 
A доери T 
ш.а а/ы; en 2) ИШ 
дх в dt T, R Te ot Т, dt 





The assumptions necessarily made in deriving equations (2.1), 
(2,2), and (2.3) are: 

(1) Zero thermal conductivity in the direction of flow Гог 
both the gas and the matrix and infinite thermal con- 
ductivity in the matrix transverse to the flow direction, 

(2) The pressure drop in the device is neglected, 

(3) The gas obeys the perfect gas lav, 

(h) The gas transport properties and the matrix properties 
are constant, 

(5) The cross sectional area of the device open to flow, Ae 
is constant, 

In addition to these assumptions it is necessary to assume that 
the heat transfer area per unit volume of the pulse tube is fairly 
large, in order that the pulse tube can be considered as a regenera- 
tore 

Boundary conditions for equations (2.1), (2.2), and (2.3) for 


the regenerator are: 


at x = 0: m = mo(t) 
lo = Tgo(t) (2.4) 
at x “< L: 1, = Tay (t) о 


Boundary conditions for equations (2.1), (2.2), and (2.3) for 


the pulse tube are: 


at x 30: m=n,(t) 
Te Too\t) (2.5) 
i x = L: 1, - Toc C6) 


Qu 





Note that these boundary conditions are identical and have been 
presented separately only to show that the position x = L for the 
regenerator is the warm end of the regenerator, while the position 
x = L for the pulse tube is the cold end of the pulse tube, when the 
Pulse Tube is operated as a refrigerator. | 

For a regenerator (and pulse tube in this case) operating in the 
Steady, cyclic manner described here, certain cyclic conditions must be 


met, At a given position x, these conditions are as follows: 


$ nAs(T, = ава Hos ol dt = 0 (290) 
òt 

Q mat = 0 | (2.7) 

g2 dt = 0 (2.8) 

О (пт + Ue T (2.9) 


Equation (2.6) says that over a complete cycle there can be no 
net heat stored by the matrix. Equation (2.7) says that there is no 
net mass stored over a complete cycle, Equation (2.8) states that the 
pressure is cyclic, Equation (2.9) results from equations (2.5) and 
(2.8) applied to a cyclic integral of equation (2.1). C, is constant 
with both time and position. The cyclic conditions, equations (2.6) 
to (2.9), are general in nature and must be satisfied by any steady, 
cyclic solution to equations (2.1), (2.2), and (2.3). 

The complete cycle will now be broken into two half-cycles of 
equal time duration. Also, the pressure is assumed to be a sawtooth 


function of time such that 





= a constant. Other approximations and 
dt, 








assumptions are as follows: 
(1) The temperature difference between the gas and the matrix 
at any point within the regenerator (or pulse Ше і5 
very small, This enables one to substitute T for T. in the 
final continuity relations, 
(2) The product of the average equals the average of the product, 
This approximation is useful in applying the cyclic condi- 


tions, equations (2.6) through (2.9). 


(3) For the regenerator only, 
9. [= СР ы 
дъ Т, T dt 
lith these conditions imposed, the final set of equations to be 


solved for the regenerator are: 


-cp = (mT ес) + ЋА Too ~ T) + cvAo 





a O 











К dt 
cd (пт_ ) + ЋА(Т__ - T) to |ЯР| то пи даи 
ра ге шие = ee 
R dt 
dm Ag |dp | 
ا‎ 2.12 
dx RT |а у == 
m(To 三 Toe) - Oy 9 (2.13) 
Е Ет (2.14) 
The boundary conditions for these equations are: 

at x = 0: те ТА 

m = m. (2.15) 


at x = L° 


" 
= 
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The final set of equations to be solved for the pulse tube are: 

















š d (.. + + С др 
С m - 
P (mT) MA, (To = 7) A - os (2.16) 
d (y rues | ple 
ер Те) + ВАТ - Т) хо ро 6200 
R dt 
dam A d Е 
о | 1-20 2 [ш | (2.18) 
ПН КА | dt 4-1) 4% Тіп KA pL Г 
n(T,. - ты) = Cy , (2210) 
Toc + Тое ео Т e (2.20) 
The boundary conditions for these equations are: 
tx = O: = 
at x * O: ТЕТ. 
mos Mo 
врх Ls Т < Те (2521) 


Note that the set of equations for the regenerator is identical 
with the set of equations for the pulse tube with the exception that 
the continuity equations (2.12) and (2.18) are different. The differ- 
ence between these two equations is that, in the case of the regenera- 
tor, the second term of equation (2.3) can be shown to be negligibly 
small in comparison to the first term, whereas in the case of the 
pulse tube, the second term is not small for all conditions of 
operation, Rea assumed that the continuity equation as written for 
the regenerator applied equally to the pulse tube, The details of 


the relative orders of magnitude of the terms of equation (2.3) for 


-11- 





the regenerator and pulse tube are presented in Appendix A. 

Equations (2.10) through (2.14) are five equations in five un- 
knowns for the regenerator and can be reduced to two equations in two 
unlmowns by subtracting (2.11) from (2.10), and then by substituting 
into this result equations (2.12), (2.13), and (2.14), plus the com- 


monly known relationship R = ср - су» The result of this manipulation 











TS: 

• ат ПА А др 

E Ti 422 == o | (2.22) 
сыл cp t 

and, 

dm А. [ар 

ор E 5 (2.23) 

dx RT |а 








The three boundary conditions (2.15) are required since, although С, 
in equation (2.22) is a constant, it is not known at this point, 
Likewise, equations (2.16) through (2.20) are five equations in 
ae unknowns for the pulse tube and in a similar fashion to that 
described above can be reduced to two equations in two unknownse 


These results are: 


‚ат һб | AQ | dp y do 1/5) [v 
二 Eois 
p “p 5 1 


4-1 (2.21) 
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and, 


A 


dm _ dp 
ах 


7 ست 


dt 


O 


Pu B 


yl 


| Y Bud 
RT 


p 5 
At Tim /| KApL 


Again, the three boundary conditions (2.21) are required, since Cy is 











unknown at this point. 


Set-up of the Differential Eguations for Computer Solution 


The following dimensionless variables are now introduced: 


MT. 
nme m/m, 
(2.26) 
= V dp x . 
To oR |а L 








Also 3% 35 assumed that the heat transfer coefficient can be 


correlated as: 


I> kmn (2.27) 


for the regenerator, and as: 
h E К (2520) 


for the pulse tube, which is assumed to be operating in laminar flow, 
When the new variables (2.26), along with the assumption (2.27), 
are introduced into the regenerator equations (2.22) and (2.23), the 
final regenerator equations result in: 
с - В а-л 
ад 2 


=a = <a 5 (2.29) 
ди monn 


2185 





dau 1 ` 
d — == 2.20 
у“ aw Ө т | | 
where 2 у= Ë os E | (2.31) 
Y | v / en, 2 
dt 





The boundary conditions for ihe regenerator become: 


at w = O: m= 1 











su. 
y |% 
at w = Wmax - ix е —) - Ey - EE е (2,32) 
T MoE To 


In a similar fashion, when the new variables (2.26), along with 
the assumption (2.28), are introduced into the pulse tube equations 


(2.24) апа (2.25), the final pulse tube equations result in: 


| + 5. (к 

СО eile хай : (2.33) 
dw "o 

спа ШЕТ |1 500 : (2.34) 
с е On 


where C is defined by equation (2.31) and 
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Jeo 
ler ии 
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The boundary conditions for the pulse tube become: 





at w “ O: m= 1 

O=1 
EST Told (2.36) 
at WS Woo, TAR a Т. ce •3 


Method of SOLUTION 


The method of solution for the regenerator equations (222097 2а 
(2.30) is to treat C as a parameter. “his makes the problem an initial 
velue problem, since both C) and m are known at м = 0. To obtain a set 
of solutions, one must know the ratio of the specific heats; $ , and 
the type of matrix used in order that n can be assigned a value. The 
parameter C is then assigned a range of values, and the equations are 


numerically integrated from w = 0 to w= Solutions for n = .59 


Wax e 
and 0 = 1.67 are presented in graphical form in Appendix С. 
The method of solution for the pulse tube equations (2.33) and 
(2.3) is similar to that for the regenerator. Now, however, there is 
a second parameter, (p , involved. Іп this case (b is assigned a value, 
© = Ф}; and then C is assigned a range of values for this Oz Again 
the ecuations are mmerically integrated from w = O to w = мрах• The 
result is a family of curves of © versus w and m versus w, with each 
curve having a specific value of C and each graph representing one 
value of the parameter Ф ə Solutions for various values of Ф and 
$ 2 1,67 are presented in graphical form in Appendix C. 
To use the curves in Appendix C one must Know Wray ane Өт, 


for the case of the regenerator. In the case of the pulse tuve, 


И хаха Ө са апа Ф must be known. Then entering the appropriate 


E 





graph with w = Wy, and e Өү» ог © = e. a value of C can be 
read. The axial temperature distribution is determined by following 
that value of C down to w = O, Also, once C is know, the mass flow 
into the top of the regenerator (or the bottom of the pulse tube, as 
the case may be) is determined by entering the graph of m versus w at 
МОБ пиву and reading off the value of m corresponding to the inter- 
section of the C curve and the w = Wax lines As in the case of the 
temperature distribution, the axial mass flow distribution is deter- 
mined by following the same C curve from w = Wace oo ae Co minis 
yields the mass flow distribution from the warm end to the cold end 
of the regenerator or from the cold end to the warm end of the pulse 
tute. 

From this point on equations (2.29), (2.30), (2.33), апа (2.31), 
plus the boundary conditions equations (2.32) and (2.36), will be 


referred to as Reats Modified Theory for the Pulse Tube. 
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CHAPTER ICE 


ANALYSIS OF THE PULSE TURE 


The System 
The Pulse Tube consists of a regenerator in series with a pulse 
tube, as shown schematically in Figure III. This, then, is the system 


which will be analyzed, 


FIGURE ITI 


WORKING SCHEMATIC OF THE PULSE TUBE 
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Analysis 
If one considers the control volume around the entire Pulse Tube 
as shown in Figure III and writes an energy balance for one cycle of 


length T, the result is: 


mol oe 











Т сер At, * Opener | у ee ИР Tep 


W 
E Q T = О , cE 


where the subscript W refers to the warm end of the regenerator. In 
equation (3.1) the first term represents the enthalpy carried by the 
gas into the control volume at the warm end of the regenerator,. the 
Second term is a known heat quantity introduced into the control volume, 
the third term is the heat associated with conduction and radiation that 
is absorbed by the control volume, the fourth term represents the enthalpy 
carried by the gas out of the control volume at the warm end of the regen- 
erator, and the fifth term is the heat removed from the volume Vye 

Now, if the control volume around the pulse tube only is considered, 


one can write an energy balance for one cycle of length T as follows: 


° ° " 
Т лы), +Q E Tp EN] 7 Q T =0, (3.2) 


where the subscript C refers to the cold end of the pulse tube and 
cae are the radiation and conduction losses associated with the 
pulse tube, 

If one now makes the assumption that the period of compression, 


J 15 equal to the period of expansion, T. and, furthermore, that 


either of these periods is equal to WT, where wis the percentage of 


=l 0a 





the total time that is devoted to compression or expension, then equa- 


tions (3.1) and (3.2) can be written as: 


Че 


ЧИ Cp (re u T W i Cheater ü losses E QH -0 و‎ (3412) 


and 


* 1 Е 
г, š tee | А + 0 уоѕѕеѕ “ дн 7 9 ° (3.2а) 


ще 


Y Cp 


By combining equations (3.1а) and (3.2a) one obtains the result: 


Ор = Yo, | mc г, - P 275% e - teh ; (3.3) 


where 


Q 


6 : 
R heater 5 а ? 


(3.1) 


and mc. | Тос = ne) is the difference in enthalpy flux at some 
point in the regenerator or pulse tube. This enthalpy flux can be pre- 
dicted by Rea!s Theory as presented in Chapter II. 

The difference in enthalpy flux at the cold end of the pulse tube 


and at the warm end of the regenerator will now be found, 


From equation (2.31), 








още ос = | E ; (3.5) 
1 4-1) (КАФ а% 

and from equation (2.13), 
сү=т т E Tze] ° (3.6) 


ЕСЕ 





Now consider the regenerator. At the bottom, or cold end, x = O апа 
m= m. = Me At the top, or the warm end, x = L, Tš Ty, and Rs m. 


Fron equation (3.6), 


W E = al: al, (3.7) 


Then from equation (3.5), 
C E БЕЛЕ _У | г 
В 4-1 КА = 


Now consider the pulse tube. At the connection between the cooler 


dp 


dt 


9 (3.8) 








my = my. At the bottom, 


section V,, and the pulse tube, x = 0 and m 


m = mo, and T 


or the cold end, x = L, m= T, = Te, From equation (3:005 


e 


me [Tye = бе | „= CENE CUM (3.9) 


Then from equation (3.5), 








C = у рү п | д 
a). „| | | чар. (3.10) 


Т pt 


By substituting equations (3.7) and (3.9) into equation (3.3) one 
obtains: 


Е ја | . (3.11) 


d = Yo, «| | 


pt 
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Now by substituting equations (3.8) and (3,10) into equation (3.11) the 


result is: 


Dr = 2 Y Cp ЕЗ 


а о 
y E "| 


dt 





=] ы-і (аш) o 
FONS 
ee a (3.12) 


Equation (3.12) is identical with equation (3.l) of Rea (lL) with the 
exception of the W factor, which apparently was omitted by Rea. The 
method of determining this factor will be commented on in the following 
chapter. 

Equation (3.12) will now be simplified. One is able to write for 


a perfect gas: 
= 
un s (3.13) 


where m is the mass, ЈЕ this equation is differentiated with respect 


to time for a particular volume, ihe result is the general relation- 


Ship: 
y wo» d | (3.11) 
R 1, dt Т, а% 


Now consider the particular volume to be Vr, the volume of the cooler 
at the end of the pulse tube. The assumption is now made that there 


is perfect heat transfer in V Tnis means that the gas in the volume 


H° 
Уң is at a constant temperature le With this assumption equation 


(3.111) becomes: 


хола 





RT, | 4% 


From equation (3:15) it fellons that 





cr > E (3.16) 
ne RT, dt 





By substituting equations (3.15) and (3.16) into equation (3.12) and by 


rearranging, equation (3.12) can be written: 











ба = ря | = | x fa]. e Б ко | io) 
Where 
n 2% wl le Be | (3.18) 
™ Ene (КВ Ја 
and 
= IT ° (3.19) 


Equation (3617) is the governing equation for the Pulse Tube. 
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CHAPTER IV 
EXPERIMENTAL PROCEDURE 


Description of Apparatus 


Rea!s Modified Theory for the Pulse Tube was tested using an 
apparatus constructed as shown schematically in Figure IV. Tne regen- 
erator consisted of a 12"-long, stainless-steel tube, 1" O.D. by 
0.020"~-wall-thickness, packed with 0.050"_diameter lead shot, The 
pulse tube consisted of a 12"-long, stainless-steel tube, 5/8" 0.D. by 
0,016"~-wallethicknesse The bottom of the regenerator was connected to 
the bottom of the pulse tube by а 1/1" 0.0. by 0.020"-wall~thickness 
stainless-steel tube.  Fine-mesh-bronze screens were installed at the 
bottom of the regenerator and pulse tube to smooth the flow, The heat 
exchanger at the top of the pulse tube was constructed from thin-wall 
l/8'-copper-nickel tubing. The rate of flow of cooling water around 
these tubes was 0.5 gal/min. A heating coil of high resistance wire 
was wound around the connecting tubing between the regenerator and 
pulse tube as a means of introducing a known heat load into the device. 
Copper-constantan thermocouples were attached along the length of the 
regenerator and pulse tube. Pressure transducers were installed at 
the warm end of the regenerator and at the water-cooler end of the 
pulse tube, The entire Pulse Tube was suspended in a vacuum-tight 
chamber for thermal insulation, During data runs this vacuum was 
maintained at approximately 2 x 10-5 mme Hg. Helium was used as the 
working gas for the Pulse Tube. A four-speed motor allowed a choice 


of cycling rates corresponding %о 22.4, ЦЦ.7, 82.3, апа 158,0 
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SCHEMATIC OF TEST APPARATUS 


NEEDLE VALVES 
SOLENOID VALVES 





Xx << >. 
RR 
50606000 
5250605 





Ў 
< | 
VACUUM JACKET | 


pulse tube 


“~— SCREENS 


PRESSURE TRANSDUCER 
TO OSCILLOSCOPE 


— 


On A 25 25 205 کک کل‎ 


CCC 


Кр? ЧУ ЧУ ДУ УУ У 


VARIAC 


zo 


`. 


WATTMETER 





reversals per minute, A complete description of the apparatus is 
contained in Appendix D, 
The important characteristics of the device were: 


(2) V 3.39 in? 


ръ 


(3) V4 = 0,673 113 


(4) Vary = 5.01 

(5) Мр Тр 2 0.965 

(6) €4 = 0.388 

(7) е = 10 

(8) 421.67 

(9) (KB). = 595 k BIU/in-hr-ft^-9F 


(10) (КВ), = 1.1) x 10° (m. ) «92 BIU/in-hreft°=°F 


(к В), assuming fully developed laminar flow, was obtained from the 


expression: 


- 174 ШК 
(KP) t = Са? 


where Nu 2 1.36 and В 5 3126 


V 
(КР) was obtained as follows: for a spherical matrix, pa = 80-62 و‎ 


where d is the diameter of the matrix spheres. KS was determined from 


Rea's experimental results (lh). 


Eos 





Method of Testing 


Basically the objective of the experimental work was to obtain 
data from a Pulse Tube that could be compared with the refrigeration 
effect which Rea's Modified Theory predicts for that Pulse Tube, 

One of the most important experimental results was the axial 
temperature distribution, particularly of the pulse tube. From the 


equation (2.26), w is defined as: 





w š V dp = e | (hel) 
L 


T¿M¿R dt 


For the pulse tube, from equation (3.15): 








° • V dp 
m = H Nm — (3.2) 
TyR | dt 
Therefore, for the pulse tube: 
V 
porre pu ees (h.3) 


Die , 
Уң Lot 

where x is the distance along the pulse tube measured from the top as 
shown in Figure III. For a given pulse tube, equation (L.3) allows 
one to specify locations along iis length as a function of Wot? This 
was helpful in that thermocouples could be located at any distance x 
along the pulse tube and that position could be translated easily into 
a corresponding wpt。 

The assumption was made that the temperature of the gas inside the 
cooler was uniform and at the same value as that of the cooling water. 


In reality this is not possible, but by using thin=wall cooler tubes 


and by allowing a short dwell time between the compression and expansion 


“26. 





parts of the cycle, "perfect" heat transfer could be approached. Tn 
this case (Т0) 4 became simply Ty. Then by measuring the temperature 
of the pulse tube at each thermocouple location and by dividing each of 
these temperatures by Тү to convert to the non-dimensional temperature 
variable CÓ, the experimental temperature distribution of the pulse 
tube could be determined readily. When this experimental distribution 


was compared with theoretical temperature distributions, one obtained 


typically: 


we max 

The shape of all the experimental curves was similar to the shape of 
the theoretical curves except in the vicinity of Ways indicating that 
the laminar flow assumption for the pulse tube was fairly accurate, 

In normal application the theoretical curves would be entered knowing 

a value of Ша чаша = However, in this case, since the entire pulse 
tube temperature distribution was known, it was more appropriate to use 
(Өг) Ar as an entering argument, This gave a more accurate value of 
E than if Ө, had been used, 


The selection of the method of comparing the calculated results 


опа 





with the experimental results will now be discussed, From equation 
(3.4) it should be noted that On includes the input power from the 
heater plus the conduction and radiation losses in the — and 
pulse tube. Since both conduction and radiation losses depend upon the 
axial temperature distribution, these losses would be approximately con- 
stant for any operating conditions of the device which resulted in the 
same temperature distribution. This suggested a possible procedure for 
testing wherein a change in the measured value of oe +$Q70sses) 


could be compared with a change in the calculated value of 0» for simi- 


2 


lar temperature distributions and for identical changes in - 








dp 
` dt 
The Ab M OÉ course, had to be determined. As a first approxima- 
tion they were analytically obtained from equation (3.17) for no 
heater input. This approximation assumed that the analytical express- 
ion is correct. This assumption may not be exact, but should suffice 
as a first approximation. The comparisons suggested above could be 


2 


made on a plot of On versus <р as shown below: 








Фу 8, 


5. 
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A smooth curve cannot be legitimately drawn through the measured points 


of the Qs versus 





al plot, since there is a third variable » Ф , 
which was not fixed, Another method of conducting the experiment would 
have been to keep ф constant, and then compare changes in experimental 
Qa with changes in predicted Qs. 

Pased on the above discussion an experimental procedure was 
selected as follows: First, a cycling speed was arbitrarily selected 


from the four speeds available, Second, a Oc was arbitrarily chosen. 


The selection of Ө, had one limitation, that being the value selected 


had to be greater than the lowest possible ©, obtainable with the 


С 
device under the limiting conditions of operation selected, This limi- 
tation was necessary so that a range of operation was available wherein 


ihe same E. could be maintained with an increasing heat load, Third, 


a ratio of p was selected which would maintain the chosen Ө. 


паж/ Ризи 
Imitially there was no heater input, Under these conditions the 
refrigeration effect consisted of conduction and radiation losses only, 


Yt took 60-90 minutes for the temperature distribution to become steadyo 





SP required to maintain the 


Once steady state was reached, the 3t 





selected =F was noted. 








Next, the heater power and were increased in such a manner 


de 
t 
that the O. remained constant. It was found that response time to 
slight changes in pressures was considerably less than the response 


time to changes in the heater input. Therefore, for rapid stabiliza- 


tion, the heater power was increased a fixed amount and then the E: 





was adjusted to achieve the desired Ө,. Once heater power was intro- 
duced, it was impossible to maintain the identical regenerator tempera- 


ture distribution, although Ө, did remain nearly constant, This, of 





course, infers that tne Qlosseg were not exactly constant for a series 
of data runs, The temperature distribution of the pulse tube on the 
other hand could be very closely duplicated. Once the steady condition 


was achieved, the P and Qheater Were recorded, By adding the 





ees to the A the Qn was obtained. By continuing this process 


а series of DEP Rudi aee points of Qn versus E4 2 could elfobtained? 
In the above procedure 2 was determined by measuring 
wattage and presented no difficulty, Because of the strong dependency 
dp 
dt 


of the Pulse Tube performance on ; provisions were made to measure 








the pressure at both the warm end of the regenerator and water-cooler 
end of ihe Pulse Tube. To check the pressure drop through the Pulse 
Tube, the water-cooler end pressure transducer was first calibrated 
under static conditions using previously calibrated pressure gauges 
which were installed at the warm end of the regenerator as a reference. 
The pressure was then allowed to cycle. Photographs of the oscillo- 
Scope traces of the pressure tranducer output were taken for different 
‘pressure settings on the gauges. From these photographs and the 
pressure transducer calibration curve, the maximum and minimum pressures 
at the water-cooler end of the pulse tube were obtained. For most 
pressure settings the measurements taken at both ends of the Pulse Tube 
were essentially identical, For cycling speeds up to and including 
82.3 reversals per minute, the only measurable pressure drop for the 
data taken was approximately 2 psi at a supply pressure of 150 psige 
For these runs it made no difference where in the Pulse Tube the 
pressure measurements were taken. With the highest speed of 158,0 
reversals, the pressure drop ranged from 3 psi at 50 psig to 2l psi at 


112 psige For these runs the p... and P min Were measured on the pulse 


es 





tube side of the device, 

Typical pressure transducer traces are shovm in Figure V, From 
the traces it can be seen that the total period can be obtained pre- 
cisely. However, the values for 8 and | P| are not quite so 
obvious, In the analysis of the Pulse Tube it was assumed that 
m = T. = WT. In keeping with this assumption the following values 


of WT were decided upon by taking an arithmetic mean of T and 7 ; 
c e 





Reversals Min WT (sec) ЧИ 
22 olt 1.0 0.372 
41.7 0.5 0.372 
82.3 0.27 0.372 
158.0 0,11, 0.372 


Because of the difference in slopes of the compression and expan= 
sion parts of the cycle and the slight non-linearity of the exhaust 


history, average values were also used to determine ЕЗ With a simi- 





lar pressure history, Rea (l) reported good results by taking: 





др ED c 
% ag = , (15) 
Atavg 
vhere 
At vg EUT = i іт. + т.) (1.6) 
апа Виа - maximum cycle pressure 
ET = minimum cycle pressure. 


A similar relationship for = was used in this case, 
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FIGURE V 
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То determine a predicted value for Qn» equation (3,17) was used: 








o 2 2 
Ор = 2 y NEG | 
| ЕЗ КЕ (ка), 
Н Plot 


where of is defined by equation (3.18). For any Pulse Tube € зи En 


and V; would be known, Y and T. would also be known depending on the 


др 


gas used and the temperature of the cooling water, МИ and q | were 








determined by the methods previously discussed, (K P Ја was known 
once the pulse tube temperature distribution became known. To deter- 
mine a value of C pt a value of the parameter ф had to be determined, 


From equation (2.35): 


> 


In the above equation Y and | ка. | were known; T, corresponds to Ty 
and was also known; At was taken as YT, which was obtained from 


equation (1.6); p was taken as Payg Which was defined as: 


+ Pmax ^ Pmin 
2 


: (4.7) 


Pavg 7 Pmin 


Thus 中 was obtained, Now since ф sul га Land (Cc) erp Were all 


pt 
known, the appropriate plot of © versus w for the pulse tube was 


entered and a Cot extracted. This in turn led to a value of m/my 


by the procedure described in Chapter II. From пути» Me was found, 
This led to a value for (KB )pe The next step was to determine Пе ee 
which was obtained from: 


Еда 





E. v 

(м). 5 КОБЕ ОНЕ : (4.8) 
s ща 
mH 


Equation (1.8) follows from the definition of w,,,, equation (2.26), 
assuming that n is constant along the connecting tubing. All quantities 
in equation (1.8), then, were known. Now knowing е апа ах)» а 
value of Ср ма5 extracted from the regenerator curves of Q versus w. 
From equation (3.18) & was calculated, since all necessary quantities 


were known. Finally, o was calculated from equation (3.17). 


Sample Calculation 


Suppose one wants to predict the refrigeration for a Pulse Tube 


with the following specifications: 


cooling water temperature = 502°R 


Ty 


Ty = 477°R 
Трън = 0,673 im 
Vp 2 3.27 in? 
Refrigeration is desired at 37678. Hence, To/Ty = 0075 апі Т/Т = 1.27. 


Helium is to be used as the working gas. The cycle rate is to be 82.3 


reversals per minute. The pressures to be used are: 


Bc 60 psig 


Patm = 15 psia 


ie 





Ф is calculated from equation (2.35) to be 2.116. From Figures 


XVIII through XXVIII in Appendix C with Mn 5.0 and Tc/TH E 0.75, а 
cross curve of 中 versus Cot can be obtained, Entering this cross 


curve with Ф = 2.16, a value of Cot = 0.117 is read. In an analo- 


gous manner a cross curve of m versus Ф can be obtained, From this 


curve mc/my = 5,8 is read, Then, from equation (4.8), 


(лак) = lel • 


Now from Figure XVI of Appendix C with w = 1.1) and Ty/Tc 271407, 
Cp = 0.85 is read, Suppose also that œ has been calculated to be 
0.00811., then XC, = 0.0069, For a cycle rate of 82.3, WT = 0.27 sec 


and W = 0.372. From equation (4.5): 

др | = | 

dt 185 psi/sec е 
In addition (KB), = 11.6 BTU/in-hr-ft?-°F, 


Then, from equation (3.17): 


Ор = 1.25 watts 





CHAPTER V 


RESULTS 


Presentation of Resuits 

Experimental data was obtained using the apparatus and the pro- 
cedure described in the preceeding chapter. Three series of data runs 
were made, Tn the first series different speeds and widely variant 
back pressures were used for the various runs. For all runs in the 
second and third series a single speed and a nearly constant back 
pressure was used, Tables I, II, and III summarize the results for 
these series of data. Figures Vi through VILI illustrate the compari- 
Sons between the predicted and the actual refrigeration loads. Observed 
and theoretical pulse tube axial temperature distributions for the three 


series of runs are presented in Figures IX and X, Figure XI shows a 


plot of m and Cot versus ф for w = (илах рё а= (Әс) еу. 
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Discussion of Results 

From Tables I, II, and III it should be noted that “XC, 15 less 
than 10% of Cot in all cases for this data and for practical purposes 
can be neglected when calculatirg on The data indicates thst es long 
as the regenerator used in the device is an efficient one, its physical 
dimensions relative to the pulse tube dimensions can be smaller than 
those of the regenerator used in this experiment before its effect would 


have to be considered. If the <C, term of equation (3.17) is neglected, 


that equation becomes: 














e 2 ` 12 | 
dg = 2y а dp ҮН 5% m 
or, for a given pulse tube, 
e ie 2 
и m (5.2) 
: 2 
where Sot is a function of ф . АП the Figures of On versus | 92 m 








this chapter are plotted using equation (5.2). As a result, all lines 


dp 2 


of ф = constant on the Qn versus 


Plots are straight lines pass- 








ing through the origin, 
For a particular series of runs a single temperature distribution 


was used in determining a C_, for any specific run in that series, 


pt 
since all runs in each series had practically identical pulse tube 
axial temperature distributions. The observed and theoretical pulse 
tube temperature distributions are shown in Figures IX and X, The 
theoretical results for each series are plotted as single curves, 

lc 


since the deviation between runs in any one series was less than 5%» 


It is also interesting to note that the experimental curves are in good 


Eo 





agreement with the theoretical curves, indicating that the laminar flow 
assumption for the pulse tube is a good one. 

From Figures VI through VIII one can determine approximately the 
accuracy of Rea!s Modified Theory. However, the accuracy shown in these 
Figures is no better than the accuracy of the variables involved in 
arriving at these results. From equation (5.2) it can be seen that бы 


ар 2 
Go 


depends on 








and C поље The quantity С + is, in turn, dependent 


upon ф e dn determining = and Ф gy average values of pressure and 








time vere used, The methods of determining these average values seem 
reasonable. However, other methods which would have resulted in differ- 
ent average values could also be justified, The importance of these 
results, then, lies not in attempting to show their accuracy, but in 
showing that the trend of the results can be explained by the theory. 
The dimensionless parameters which evolve from the theory are most 
veluable in understanding the problem. For example, consider Runs 2 and 


5 from Series 1 of Figure VI. Note that the Е: : is approximately the 
Е 


same for these two runs and yet esper for Run 2 jis more than twice 


Q. for Run 5. This can be explained, of course, by the fact that 
exper 

for Run 5 the value of © was almost twice as large as (p for Run 2. 
AS a result Cot for Run 2 was more than one and one-half times as large 
as Cot for Run 5, Ап indication of how С 


is shown in Figure XI. hus, Qs, for Run 5 was approximately one- 
cor 


varies with ф in this case 


да A 2 e di 1 For 
helf as large as рео for Run 2. The large difference in ф Го 
these two runs is attributable to the fact that the back pressure for 
Run 5 was four times larger than the back pressure for Run 2; also 

Run 5 had a WT of 0.27, while YT for Run 2 was 0.50. The effect of 


these quantities on Ф cen be easily seen by examining equation (235) 


Е 








Another example of interest are Runs 5 and l| of Series 1. Note 
that the total experimental refrigeration for these two runs is about 
equal, yet ЕЗ! is more than twice as largo for Run 5 as for Run Ше 
The explanation, of course, was that @ for Run 5 was almost twice as 
large as Ф for Run h, The values of 中 were different, since WT for 
Run 5 was 0.27 opposed to 0.50 for Run l and the back pressure for Run 5 
was 20% greater for Run 5 than for Вип h, This resulted in Cot for Run li 
being more than one and one-half times as large as Cot for Run 5. This 
difference in Cot due to Ф наз enough to counteract the opposite effect 
on баео due to 2 and, hence, both conditions resulted in pro- 


ducing approximately the same refrigeration effect, 


2188 





CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 
Conclusions 
From the analytical model of Chapter II, there appeared three 


principal parameters which characterize the operation of a Pulse Tube, 


These parameters are: 





hcl ше a l ‚ (2.31)‏ م 
г | 7 ел, |%)‏ 
dt‏ 

W = у ар 9 (2,32) 


which were described by Rea (Nh), plus a third parameter Ф » defined as 


ЫК , (2.35) 


KA pl At D 
The parameter Ф resulted from the fact that the ratio 


1 зр //h 27, 
— LL — — | for the pulse tube was much closer to unity than 
T. òtl/ ITE òt 
g 8 
for the repgenerator, Thus, the simplifying &8sumption which neglected 
the term [P 21, in the continuity equation (2.3) in the case of 
ще д +, 


the regenerator could not be made in the case of the pulse tube. 


From the theoretical and the experimental results obtained from 


thé Pulse Tube under consideration, i 
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Theory as presented in Chapter II does describe the performance of a 
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typical Pulse Tube, There was some deviation between the refrigeration. 
effect measured and the refrigeration effect predicted. However, this 
deviation falls vithin acceptable limits for heat transfer data in the 
cyclic case, | 

It was apparent from the data obtained that there was good agree- 
ment between the experimental and the theoretical pulse tube axial temp- 
erature distributions. This leads one to the conclusion that laminar 
flow in the pulse tube is a good assumption. 

From the analytical results obtained from the Pulse Tube, it was 
concluded that the effect of the regenerator on the overall refrigera- 
tion was negligible. Consequently, the regenerator could have been 
made smaller in this case before its effect would have become signifi- 
cant. 

Experimental and theoretical results indicated that there is a 
unique temperature distribution for a pulse tube operating with a 


specified Oc. 


Recommendations 

There are several areas of Pulse Tube refrigeration which should 
be examined. First, the uniqueness of the axial temperature distribu- 
tions of several pulse tubes with different values of Her operating 
with a specific Өс should be investigated over a wide range of ф ° 
If this uniqueness can be demonstrated for a significant range of w,., 
and ф , then a suggested method for presenting the computer solutions 
for the pulse tube differential equations is as follows: 

C0) Flot for a fixed W ay ê mr of cross curves of Cot versus 


ф for various values of = 


(2) Introduce theoretical axial temperature distributions for 


Ва 





some value of Ф which covers the maximum feasible range in w. 
With the temperature distribution curves described above, the е end 


Bex of any pulse tube can be used to determine a curve of constant C t" 


By following along this constant C ¿ Curve to the value of w for which 


p max 


the cross curves of Cy. versus Q apply, a corresponding value of е 
can be obtained. This value of Oc represents a single curve on the plot 
of the cross curves of Oy. and Q. Thus, the relationship between Cp. 
and 中 becomes known, 

Second, a method for determining an optimum size regenerator to be 
used with a given pulse tube should be devised, 

Third, an optimizing procedure for the design of a pulse tube using 


Rea!s Modified Theory should be developed. 
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APPENDIX A 


REAIS 2-PART MODEL AND TEE DEVELOPMENT Or 


Eee eee ee eee a TT rt 


IHE GOVERNING DIFFERENTIAL EQUATIONS 


General 

As a mathematical model for his analysis of thermal regenerators 
subjected to rapid pressure and flow cycling, Rea devised his 2-Part 
Model. From this model Rea developed the differential equations that 
apply to regenerators subjected to the above conditions, Rea proceeded 
to apply certain assumptions which simplified these equations. The 
final result was a system of two differential equations, which Rea 
solved numerically. He then reasoned that since a pulse tube is 
Simply an inefficient regenerator, the same differential equations 
should apply to it. However, not all of the assumptions which led to 
the regenerator differential equations hold true in the case of the 
pulse tube. Thus, the final differential equations which govern the 
pulse tube are not the same as the differential equations which govern 
the regenerator. This Appendix is broken into three sections. Ја 
Section A.1 Rea!s 2-Part Model and his development of the differential 
equations for the regenerator will be presented in detail as in (li). 
Section A.2 deals with the development of the differential equations 
for the pulse tube. The latter section will not be presented in as 
rigorous a manner as Rea!s development in Section A.l, since the two 
developments are similar. All assumptions made in each case will be 
described clearly. Section A.3 discusses the method of solution for 


the regenerator and pulse tube differential equations. 


Zen 





Кој 


REA!S 2-PART MODEL AND DEVELOPMENT 


—— | Ci E E REI B A A i ee 


Rea!s 2-Part Model, which is shown in Figure II of Chapter I, is 
defined as a regenerator in which the mass flow rates at each end of 
tne regenerator are in phase, By "in phase" it is meant that during 
each half-cycle the flow is always unidirectional within the regenera- 
tor, i.e. during the pressurization half-cycle, mass flows in at the 
top of regenerator and exits at the bottom, and during the expansion 
half-cycle, mass flows in at the bottom and exits at the top. 

The operation of Rea!s 2-Part Model is as follows: the inlet 
valve ovens admitting high pressure gas into the regenerator. The 


9 


pressure increases within the regenerator and the mass flow rate, nos 


exits at the bottom. The mass flow rate at the top, - although 
unknown, is certainly greater than m. since some gas must enter and 
occupy the volume within the regenerator. When some maximum pressure, 
Pax? ÍS reached, the inlet valve closes and the exit valve opens. 


During the expansion half-cycle mass flow equal to m, enters at the 


O 
bottom and mess flow equal to mn exits at the top, During this half- 
cycle the pressure decreases to its minimum value, реа Опсе а 
steady, cyclic state is reached, the above processes are repeated con- 
tinuously with the temperature at the top equal to Sy and the tempera- 
ture at the bottom equal to T ° 

t is interesting to note that, since m. is always greater than 


ә 


2 all the gas entering at the top does not reach the bottom, Thus, 


c EE 





it is possible for some gas to remain always within the regenerator 


and flow back and forth as the pressure cycles. 


FIGURE XII 


WORKING SCHEMATIC OF THE REGENERATOR 


in E DS Se шш э хш = 0 


Control 
Volume 





Development of the Governing Differential Equations for the Regenerator 


The development which follows with regard to the regenerator 
differential equations is entirely attributable to Rea (h). 


Rea!s initial assumptions made concerning the development of the 


regenerator differential equations are as follows: 
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(1) The pressure drop in the regenerator is neglected, 

(2) The gas obeys the perfect gas law, 

(3) The gas transport properties and the matrix properties are 

constant, 

(Ц) Тһе cross sectional area of the regenerator open to flow, A, 

is constant. 

Figure XII shows a schematic of the system to be analyzed. The 
system consists of a regenerator of length L, which is filled with a 
matrix material of large heat capacity, e.g. spherical lead shot. The 
subscript m is used to denote the matrix. The void volume within the 
regenerator is denoted by V, and is equal То А е During each cycle 
cf operation the gas in the regenerator is alternately compressed and 
expanded E the pressures p end D e 

max min 

In accordance with the assumptions listed previously, energy equa- 

tions will be written for the gas and for the matrix contained within 


the control volume of length dx, as shown in Figure XIl, 


Gas Energy Equation: 


Energy flux into - Energy flux out of Rate of Change of 
the Control Volume ” the Control Volume * Energy within the 
Control Volune 





па + |- код әт = па + | Ба] ax - кА д1 + КА 9 | 2 Тр dx 
| 2.24) Е E g'o ен 


* 


D Р 
+ sl T. ` T йк + 2, шы : 





or, upon simplifying the above equation; 


2 
3 • à lo u d = 


2578 





In equation (A.l), the first term is the change in enthalpy flux 
across the control volume, the second term represents the gradient of 
wne energy transport by gas conduction, the third term is the сопуес- 
tive heat transfer between the gas and the matrix, and the fourth term 


15 the rate of change of energy within the control volume. 


Matrix Energy Equation: 





merey flux into КЕ Enerey flux out of » Rete of change of 
the Control Volume ~ the Control Volume Energy within the 
Control. Volure. 


or, upon simplifying, 


> T 
ща що а + Ат lo ~ T - = Mea òT е (До?) 
x? àt 


In equation (A.2) the first term is the gradient of the energy 
vransport by conduction along the matrix, the second term is the con- 
vective heat transfer from the gas to the matrix, and the term on the 
right is the rate of change of energy of the matrix. 

Equations (А.1) and (A.2) are energy balances for the gas and 
matrix, respectively, Subscripts g and m refer to the gas and the 
ratrix, respectively; h represents enthalpy and u represents internal 
energy; hg is the heat transfer coefficient (which has not been assumed 
constant) and k is thermal conductivity; 4 is the total cross-sectional 
area of the regenerator and A, is the cross-sectional area open to 


ЗЕ 





пате»? ае) 





flow (A, ШЕ А): An is the heat transfer surface area of the regenerator 
y + 3 2 А: ° ° ° ar 
per unit length (ft^/Tt); and M is the matrix mass per unit length (15/12). 
“he gas continuity equation will now be developed, 
Gas Continuity: 
rass flow into the _ Nass flow out of the Rate of storage of 


Control Volume = Control Volure mass within the 
Control Volume 


МАТЬ Sr А 
п 2 | п +90ах| + 319 4 ах | 
| ax | 2| go j 
me, 228 : (аса) 
Ó x O òt 


Equation (A.3) is the continuity equation and it is perfectly general 
(within the assumptions listed previously), as it applies to either 


steady or unsteady flow, With the use of the perfect gas law, 


Ely 
equation (A.3) becomes 
Ó x R gt To 


Equations (A01), (4.2), and (4.3) along with appropriate boundary con- 
ditions are the governing differential equations for any type o? regen- 
erator operation. Equation (4.1) will now be simplified by using 
equations (A,h) апа (1.5). 


For a perfect gas, one may write 
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h = на T Cote (1.6) 


and 


с 
1 


BE _ о КА 2 


where ho is the gas enthalpy at some reference temperature, 
The first term in equation (A.1) becomes, with the use of equa- 


tions (4.5) and (4.6): 





ај аа à [s Br оно 9 
ax 8 у 


(4.7) 


The fourth term in equation (A.l) becomes, with the use of equa- 


tions (A4) and (4.6): 





À V, 
A, ò | съ | = ойо Ваш o к” Cyto op А (4,8) 
Rue в 2\7, R dt 


The substitution of equations (A.7) and (4,8) into equation (A.l) 


yields: 


o 2 = 
c. 93 |mT л рза от) зава -0 . (A9 
s ; SED X2 | В òt cz 
Equation (4,9) is the final form of the energy equation for the gas. 
Equations (A.9), (A.2), and (4.5), along with the appropriate 
boundary conditions, completely describe the performance of tnermal 
regenerators consistent with the assumptions listed on page 57. 


The assumption will now be made that terms involving thermal 
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conductivity are negligible in comparison to the other terms in (A.9) 
and (A.2). That this assumption is consistent with some practical 


cases 1s shown by taking some typical valuesr 


kg 2 0,07 PTU/hr-ft-°P, A, = 020019 #52, ue Dg nU 





«= kas 
An = heli3 тъ? тъ, ће 2 235 BTU/hr-ft^-Or, Oo s d аре 150°F/ft“ 
5х 12 
Then: 
e1 
KA, 228 = 0.020, BTU/hr-ft, and 
à x° 


hrAn (Te = ново ВЈ ће се 


Thus, the gas conduction term is approximately 0,001% of vhe corn- 
vective heat transfer term of equation (4,9). 
in the matrix energy equation, (A.2), with nc 20 PTU/hr-ft-"F 


ard A = 0.005 +2; 


= 
к> 

^ 

к? 


3 
o 
PS 

ho 


= 15 BTU/hrefto 


Therefore the conduction term for the matrix is less than 1% of the con- 
vective heat transfer term. Then neglecting axial conduction, the final 


basic differential equations become: 





à m | Cio д 2 di 

c, S [mT| + пина т | meteo (4,10) 

Pax | з чч 8 R òt 

haa з, - z| = Me, ó (4,11) 
8 àt 





Cy 
3 © 


| 


9 
+ 


(405) 


I 
1 
m 


C 
Ра 
o 


In equation (A.10), the first term is the gradient of the enthalpy 
flux alcng the regenerator; the second term is the energy flux due to 
convective heat transfer; the third term is the work of compression. 
Equation (A.11) means that the increase in internal energy of the 
matrix equals the energy given up by the gas by convective heat trans- 
fer. Equation (A.5) is the Еее equationo 


Boundary conditions for the above differential equations are: 


KOs m= m(t) 


EJ 
if 


g ^ Tgo(t) (4,12) 


ар ха Г: а $ 


Samolaticacion of the Differential Equations 





Equations (A.10), (A,.11), and (A,5) are three equations in three 
unknowns: 7, 1.» and n. Inu uc puc function of t is known, then 
solutions consistent with the boundary conditions, equations (4.12), 
will exist. Unfortunately, equations (A.10), (A.11), and (A.5) are a 
system of differential equations for which no general solution is known. 
The method of approach, then, will be to determine solutions which 
satisfy certain conditions over a complete cycle, but which are not 
exact in every detail at each instant within the cycle, 

A cyclic integral of equation (A.11), while holding x constant, 


yields: 


- 





г гл 
9 hin y - т пре us, or (4,13) 


The above follows from the fact that once steady, cyclic conditions 
have been established, the matrix temperature at a given point in the 
regenerator oscillates around a value of the temperature which is г 
function of x only. From this point on, h will be used for the heat 
vrensfer coefficient rather than Ne 


A cyclic integral of equation (A.10) with x constant yields: 


A 《过 | :zs | dt + бт Ep - т) dt + vedi dt = 0, (4.11) 
= © 


Buty T Go = 70) by бегаа оп, Sance,p 1S a function of t only 
о 
(no pressure drop in the regenerator). Then by using equation (4.13), 


equation (A.11) becomes: 


§ 2 Ка dt z 0, (1.15) 


опис, will De continuous in x, the order of differentiation and 
о 
integration can be interchanged and equetion (A.15) becomes: 
д (ir, | ee 0 
n - = 
дх 5 
Gis 


4 (i, dt = C. (for all x). (4,16) 


For the 2-Pert Model, the vressurization part of the cycle is 
denoted by the subscript c and the expansion part of the cycle by the 
subscript e. The mass flow rate at a given point within the regenerator 


will be cyclic, passing through zero at reversal. The mass flow versus 


~ 


ro 





time might appear as: 


A 





Since there can be no net mass storage within the regenerator once 


steady, cyclic conditions prevail, 


T > Q (ръж 


and the positive area must exactly cancel the negative area in the 

above diagrame Similarly, T will be cyclic with time, but in some 
© 

unspecified manner and т. dt will not, in general, equal zero, 


Then, the product ml, might have a shape as follows: 





° 
where ihe $ an, GD em О from equation (A.16). This means that the 
positive and negative areas in the above sketch do not, in general, 
0 
cancel each other. The assumption will now be made that mI. versus 


time can be represented by the dashed lines in the above figure, 


E 





Then: 


Ó iin at = ES < t = Co 
C e 


ar led eT EE с. s (A.18) 








where T is the blow time during compression and To is the blow time 


during expansiones 
Furthermore, if (s = E (compression period equals expansion 


period), equation (A.18) can be written: 








o o 
mer + m T 


E e ge^ C5 9 (2.19) 


where 


Cy = Ue ae = а d 


Fron equation (4.17), 


9 ІСТЕ (icat s ра 0 4 
с е 


and, | mat = E : (4,20) 
с е 


With tne assumption that the m2ss flow rate can be represented by an 


average flow rate which is constant with time over each half-cycle, 


equation (4.20) becomes: 
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and again letting ЕЕ ЕЕ 
Man, қ | (4,21) 


Equation (A,20) and, hence, equation (4,21) holds at any given x 
along the regenerator. From here on, from the result of equation m T 
nm will be used without subscript to denote flow in either the compression 
or expansion half-cycle and m is a function of x onlye 


With the assumption that 


mre (A(t), 


equation (4,19) along with (A.21) becomes 
Mm (Too = Tye) = Cy (4.22) 


Equation (4,22) is a very important and useful result. There are no 
restrictions on C4, Jt can be positive, negative, or zero, depending 


upon whether T is greater than, less than, or equal to Тез? 


gc 
The usual heat transfer correlation is assumed, 


h = K z А 


where n is a positive exponent vhose value depends on the particular 
regenerator matrix. In general 0.5 = п © 0,7. Hence ‚his either 
constant or a function of Е only. Using equation (422105 hwill be a 
function of x only end will be written as he Then from equation 


(2015): 


hA 2 сал ф (p - Da: 
пас T) dt TII T) 50 ; 


since h is independent of t, or, 
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(T, Se) ates | По азот o ° (22%) 
е 


plot of (To - T) versus t might appear as: 


N 


(T, - T) 





From equation (A,23), the positive area on the above sketch must equal 
the negative area, As before, assuming (7, ~ T) can be represented by 
Т 


a 


some average value and letting |, = е 


, equation (A,23) becomes: 


Кет 2 (Тт (А221) 


e 


oince the blow-time, f , is usually short, a reasonable assump- 
tion to make is that the temperature of the matrix at a given x will 
vary very little during a half-cycle and, in fact, T can be assumed 
TI time. In that case, TS T(x). The validity of this 
assumption can be demonstrated by the following example. From equa- 
tion (A.11): 


h Ap (Tọ - ?) = Mem 51 (4,11) 


For a regenerator 12" long and 1" diameter packed with lead shot of 
0,050" diameter, Me L z ,068 BTU/?F, Taking h z 235 BIU/hr-ft^- P, 
Ат = 1.0113 ft^. and assuming (T, - T) 29р, 

hA, (T, - T) с 2090 BITU/hr - 0.068 21 
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end ре. = 30,600 ME = 8.5 °F/sec 


v 


With a blow time of 0,1 second, AT = 0.85 °F, During a complete cycle, 
st a given x, the matrix temperature would oscillate around 2 mean tem- 
verature with an amplitude on the order of 0,85 °F, Thus, the assunpilon 
is made that T is independent of time and is a function of x only, 

By applying the above assumption to equation (A.2h), one is lead 


to the important result that: 


3 | 


T A O 


e (1.25) 


At first glance equation (A,25) seems obvious, but a careful exanina- 
tion of the preceeding discussion will show that two key assumptions 
were necessary in deriving equation (A.25). The assumptions were 
(1) the mass flow rate, m, is a function of x only, and (2) T is a 
ПОСТО of x only, 

By carrying out the time-differentiation in the continuity 


equation (A.5), one obtains: 


à n ped T 
= = Nc ge = м — o. ° (A.26) 
à X R To ot lo gt 


The ratio of the first term in the brackets of equation (4.26) to the 


second term is: 


With AT, ШАТ = 0.85 °F  , the above ratio becomes: 


„608 





m 


СЭЙ 1521 - |50 to 625, 
А1, Е 


"hus, the second term in the brackets of equation (A.26) is less than 
1% of the first term and, consequently, the second term will Бе 


neglected. Equation (4.26) can now be written as: 





А. а 
= AL a • ( 527) 
К а: 


El Es 
er e) 


09 


Total dorivatives sre used in equations (A.27) since p is e function сї 


t only and m is considered to be a function of x only. Furthermore, T 
will be replaced by T in equation (A.27) since (T, - T) is small 


(29 F to 5% F), This assumption also enables continuity to be satisfied 


9 
over a cycle, Since m has already been assumed to be equal and opposite 
for the compression and expansion half-cycles, it is reasonable to expect 


< to be the same for both half-cycles. Then the continuity equation 


is written 25: 





Ag dp 


A 2220-72 (4.28) 
dx RT Gt 


s m + š : 
Tt is now assumed that = = ^" constent, Thus, pressure is 


assumed to be a sawtooth function of time as shown in the sketch below, 


where the zero-axis is any arbitrary pressure level: 


ct 


E OE 





Letting ар = ° constant also maintains consistency with the previous 
assumptions that m and T ame functions of x only. 

It is now possible to present the final formulation of the regen- 
erator problem. From Figure XII, the positive x direction for the 
regenerator is from bottom to top, Mass flow in this direction is posi- 
tive. During pressurization the flow is assumed to be from top to 


pottom, or, in the negative x direction. The gas energy equation which 


was derived for positive m with positive = is written for thu two 





cu 
helf-cycles: 
а > m m | PA а eo 
«+ С. mT + hûn | 1 - : “> C O 22 = 0 • (4,29) 
EE "18 8 | 





D (2220) 


In equations (A.29) end (4,30), = is written as an absolute value 
with the plus or minus sign placed in front, depending voon whether 


pressurization or expansion is occurring in tne regenerator,. The cor- 


tinuity relation, equation (As28), is written: 


: (А31) 


dm _ A, lap 
dx RT | dt 


For convenience, average value bars have been omitted. Equation (A.31) 
holds for both half-cycles, since in the compression half-cycle x is 
negative, but ср is positive, whereas in the expansion half-cycle 

е 


m is positive, but “P js negative, 
dt 


The final formulation of the problem consists, then, of equations 


Eo 





(As29), (A230), and (A.31) along with the results previously attained: 


ac. - T | S o (A.22) 
= Be L 
Tec * 7ga = 2T Я (1.25) 
Boundary conditions for the above set of equations are: 
at x 2 0: ES 
m Tm 
ас xs Le ТУЕ ZI Ф (A032) 


Equations (As29), (A.30), (As3l), (As22), and (A.25) are five 


~ 
equations in five unknowns: m, T, T Tze? апа 64. However, tinere 


пе. 
are only three boundary conditions. Therefore, it is necessary to 
eliminate the gas temperatures from the above equations. 

If the differentiation indicated in eouations (A.29) and (4.30) 


is carried out there results: 





МЕ ЕНИС (52) 
at 








c mm dm IX Е бр|- O , (A.30a) 
—9 ща R dt 








152 
‚ [ат ат ал 
æ ge مڼه‎ ge T “Р == i = 
cp m | = = | To x - + hÀ p т 3 
2 а 
+ судо = ZO Е (А 50 
R а 








25ДЕ 





Now by substituting equations (4.22), (4,25), end (A.31), plus the 


relationship R 2 cy - cy, equation (4,33) reduces to: 


DES rr MISES Т 
С 20, 2 dt я T 








which involves only ihe matrix temperature, T. 


The following new variables are now introduced: 


Ә 2 т/т 








о 
m 2 п/п, (4.35) 
W = u ар 到 à 
T m R at | L 


Also the heat transfer coefficient is now assumed to be correlated by 


h = Km : (А.36) 


where K is a constant with the units BIU/hr-ft*-°r. 
With the aid of the new variables, equation (A.36), and the 


definition X = сы , equation (A.3l:) can be written: 


C 
E 355-1 5 
а cy Re (4.37) 
vhere 
„ Х-1 КАТС 
С з E ED 2 (А.38) 


Гот 


The continuity relation, equation (A.31), becomes 


um a 


1 A, 
w S (4.39) 


2723 





Boundary conditions for equations (A.37) and (4.39) аге: 





at w = 0: ПЕТ 
Ө: і 
4. E = V code e - = 
ао W = y „ = Е — => (А.Ц0) 
тах ттк ]|аї p 








Equations (A.37) and (A.39), along with the boundary conditions 
(A0l/0), are sufficient to determine m, T, and the constant C. Actually, 
since everything in the definition of C is known except С.) determining 
C is essentially determining С, • 

Thus the problem for the regenerator is reduced to one of solving 
a set of two ordinary differential equations subject to the bouncary 


conditions mentioned. Tne method of en will be discussed later 


in section A.3. 


-73- 





А 2 


THE GOVERNING DIFFERENTIAL EQUATIONS FOR 


| — Q@& m — - - —— < то = жү — G s s wm m wag c 


Euer ome d Ee ары 


A pulse tube, as long as its heat transfer area per unit total 
volume is fairly large, can be considered to be an inefficient regen- 
erator. Thus, a similar development to that leading to the final 
differential equations for the regenerator could be presented for the 
pulse tube. However, such a development would be repetitious. The 
procedure that will be followed will be to accept the assumptions 
listed on page 57 and start with equations (A.9), (A,5), and (4.2), 
which are the gas energy equation, the gas continuity equation, and 
the matrix energy equation respectively. Appropriate assumptions 
regarding the relative magnitude of the axial conductivity will then 
be applied. Next the cyclic conditions will be presented. At this 
point the relative magnitude of the terms of the continuity equation 
will be examined. The final formulation of the pulse tube problem 
will then be possible, The pulse tube differential equations will 
evolve in much the same form as equations (4.57) and (A.39). It 
will be of some assistance to present here a schematic for the 
2-Part Model for a pulse tube. This is shown in Figure XIII, Note 
that x is measured positive from the top down, which is directly 
opposite from the case of the regeneratoro 

Equations (A.9), (A.5), and (A.2) in the form in which they 


were previously developed are as follows: 


E 





C 2 (in, - КА мали DEL + vo 0P = 0 (A.9) 
E 18 FEET 


on = = ES д P (A.5) 
Ox К Тт 
FIGURE XIII 


T 
0 
x 
Control TT 7 == 
Volume Ñ ax 
| D 
e 
Prax Р min 
када T E : = Me 21 (2) 
д х2 E T 


As in the case of the regenerator, the assumption will now be rade 
that the terms involving thermal conductivity are negligible in compari- 
son with the other terms in equations (A.9) and (4.2). Some typical 


values which substantiate this assumption are as follows: 


k. Z 0,07 BIU/hr-ft-?F, A, = 0,0019 s), — (T, - |) = 2°F, 


E 





_ 2 
A, 2 0.155 ft^/ft, — ho 2 6118 ErU/e-ft2-9p, 


= — z 120°p/zt2, 
L 





Then: 


2 
KA, 2 Ig - 0.016  BTU/hr-ft, and 
à x* 


hoe (T, - T) 2 1.92. BIU/hr-ft. 


Thus, the gas conduction term is less than 1% of the convective heat 
ransfer term of equation (4.9). In the matrix energy equation (A.2), 
with k = 10 BIU/hr-ft=-°F and A z 0.000225 ft?; 


ка 21 | 
A — = 0.270 BTU/hr-ft 


yx 
the conduction term is approximately 14% of the convective term. If the 
device were operated with a (Т, - N) of 190°F, this percentage would 
increase to approximately 22%. Obviously this assumption is not as good 
as in the regenerator case, but without it the problem would become 


exceedingly complex. Therefore equations (A.9) and (A.2) become 





respectively: 

o nT + hig (тост зае др = 0 (A.10) 
дх 5 5 R òt 

and 

hon т, ~ T| = Me, = E (4,11) 





The cyclic conditions for the pulse tube are identical with those 
for the regenerator. Thus, equation (A.22) applies equally as well to 


the pulse tube as to the regenerator. Again the heat transfer 
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correlation is assumed to ђе: 


ШЕ K. mn 1 
For the pulse tube, however, laminar flow is assumed. At the moment it 
is not intuitively obvious that this is a valid assumption. However, 
experimental data indicates that laminar flow is probably the closest 
approximation that can be made for the pulse tube. Thus, п = О, апа 


the heat transfer correlation for the pulse tube becomes: 


h = constant А 


Now, proceeding with an identical argument as that used іп the case 


of the regenerator, one obtains equation (4.21): 


МИ = ON) (A.21) 


Here the assumption is made that the temperature of the pulse 
tube wall at a given value of x will vary little during a half-cycle. 
Therefore T can be assumed to be constant with time, and T = T (x) only. 
The validity of this assumption can be shown by the following example. 


From equation (А.11): 


Dir SD ESCORT (A.11) 


For a stainless steel pulse tube 12" long with O.D. = 0.625" and I.D. = 
0.591", Mc,L = 0.011 BTU/°r, and taking h = 6.18 BTU/nr-ft*-°F, 


Anl = 0,155 ft?, and assuming (T, ME 


жэ = д T | 
ВАК T)z 1.929  BIU/hr = 0.011 Sm t 
and 


д 


ur. 


- 17h.5 “F/hr. = 0.0185 “F/sec, 


о 


t 


5775 





Thus, with a blow-time of 0.1 second, АТ- 0.0059Е. During a com- 
plete cycle, at a given x, the pulse tube wall temperature would 
oscillate around a mean temperature with an amplitude of approximately 
0.005 ^F. 

Thus, the assumption that T is independent of time and is only a 
function of x, is a valid one, 

By applying this assumption to equation (A.2l), one obtains, once 


again, the result: 


ПОЕ ОТОБУ (А.25) 


By carrying out the time differentiation of the continuity equa- 


tion (A.5), one obtains: 


xm e (4,26) 


The ratio of the first term in the brackets of equation (A.26) to the 


second term is: 


1 T T 
і др Ser а 
T, gt Tg ó t р AT, 
From equation (А.22): 
2 1. = VY (1.11) 
m 
From equation (A.38): 
= See V dp 
pM г ee 15 (А242) 
1 а окш |а 








By substituting equation (A.12) into equation (A.ll), one obtains: 
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ен 0 Je : (А.43) 
j-l m pl | dt 
Thus, 
ER ooo mat 
p AT, 26 Ущ У p Ë ` 


Now using typical values of Ф = 1.67, KA qL/N = 11.6 BTU/in-hr-ft^-?F, 
ШЕ 50; m/mo = 3; At I 27 sec, z 75 psia, Tg = 390°F, the above 


ratio becomes: 


A 


p 


=3 


г; 
[> 


lg 

For the regenerator this same ratio was 30 to 625, and the second term 
was neglected. Rea!s theory also neglects this term for the pulse tube. 
However, the above example indicates that the second term cannot be 
neglected in the case of the pulse tube. Fquation (4,26), then, remains 


ass 





am A /l dp p ат 
O = N (A oli) 
ах в |7, % T = dt 


Total derivatives are used because p is assumed to be a function of 
t only, Ще (at a given x) is considered to be a function of t only, and 
m is considered to be a function of x only. Also, since (T, - T) is 
О 
small (2°F, to 5°F.), T, WiN be replaced by T 11 едча оп (4.13). Then 


the continuity equation can be written as: 


dm A {la ат 
dx R та T a 


ES 





Now the assumption is made that е can be approximated by 


dt 


АТ : 
A o Once again яр is assumed to be constant for both half-cycles. 


Now by substituting equation (A.l3) into equation (A.45) and using the 


fact that 2 =- + constant, equation (А.Ц5) becomes: 

















dm А |1 | dap хрра мо | У ||ар 

dx R |T|dt (5 x ru 

or, simplifying, 

са до | др 

=== = п-в (доц?) 

ах RT] dt 

where 

и EN ње 4 : (4,18) 
5-1) дъ Тіп/) КАФ. 


In equations (A.46) апа (А.Ц7) average value bars have been omitted for 
convenience, 
As in the case of the regenerator, the gas energy equation can be 


written for the two half-cycles: 








”cp = ы + һАт Tec -T | + Cyto an - 0 ə (А 29) 
ах R dt 
d е Ãg | dP | ~ 
cm + Mm E - | COIS 0 í (А.30) 
d 可 i n R lat 








The final formulation of the pulse tube problem consists, then, 


of equations (A.Ll7), (A029), (A.30) along with the results: 
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m т. - E oM š (A.22) 


T +7 _3 oT А 22 (4.25) 


Boundary conditions for the above set of equations are, as for the 


regenerator: 
atxz0: TT (1.49) 
m = mo 
aux = 7: TS To 


Note that the condition at x z L is different for the pulse tube than 
for the regenerator only in the fact that the cold end of the pulse 
tube is at x = L whereas the warm end of the E EAS at 
x = L. Thus the subscript C stands for "cold" in the pulse tube bound- 
ary condition, whereas the subscript W stood for "warm" in the regen- 
erator boundary conditione 

As in the case for the regenerator, equations (A.29), (4.30), 
(A,l7), (4,22), and (4,25) are five equations in five unknowns: 


m, T, T ges and Cj, Since there are only three boundary conditions, 


gc? : 
it is necessary to eliminate the gas temperatures from the above equa- 
tions, 

If equation (A.30) is subtracted from equation (A,.29) and if 
equations (A.22), (A.25), and (A.lı7) are used to eliminate the gas 


temperatures, the resulting equation involves only the temperature of 


the pulse tube vall: 


• ат ГА А. 1а 1-B 1y -1 
m — = т = Ср (2-в \У -1 O ° (А.50) 
ах . 2cpm cp dt Y -1 
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Now if B is replaced by equation (A.Ll8), equation (A.50) becomes: 














O = | ле || УУ, 
ME л, o SE 2-1 E 
dx 20m Cp |dt 。 








4-1 
(А.52) 


The variables O, m, and w, equation (A.35), plus the heat trans- 


fer correlation, h = K, are now introduced into equation (4,51), 








obtaining: 
2 46 Sell. 中 
i — ano em C I -œ Ae 2 
m Ә.Е i 5 (А.52) 
where 
ф 2 2|-- RI , (1.53) 


At is Qtave for the entire cycle, and p is P for the cycle which is 


iaken as: 


_ p... = P 
Pave = Pmin ~ p (A, 5l) 


The continuity equation (A.l7) becomes in this case: 


а Фо i 
223 С . (4,55) 


The boundary conditions for equations (A,52) and (A,5h) are: 


at w Š 0: пет (A056) 
O21 
Е Се. - 6 














Thus, equations (4.52) and (A.55) together with the boundary 


conditions (A.56) describe the problem for the pulse tube, 
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АЗ 


METHOD OF SOLUTION FOR THE DIFFERENTIAL EQUATIONS 


WIE Un Lu ee mmo | фр лыты — V- — ce n s Ф mus ee eee eS 


Regenerator 


If equation (4.37) is solved for ЗӨ , 
ди 


4-1 len 
с-- = 





ав . 
n ne - (4,57) 
m n 
Also, 
am - 1 : (4.39) 
dw e 


The simplest method of solution for the set of equations (A.57) 
and (A.39) is to solve them numerically with computer aid. The 
simplest approach is to treat C as a parameter and assign a range of 
values to it. The problem then becomes an initial value problem, since 
everything on the right side of equations (A.57) and (1,39) 26 w = 0 is 
known. The independent variable w can then be run out to sufficiently 


large values to include any desired w... 


oe then becomes an unknown 
quantity, but can be determined from a plot of © versus w for a given 
value of C. On the other hand, for a given value of O,, and Ways the 
corresponding value of C can be read from the same plote 

Figure XIV shows a hypothetical set of solutions corresponding . 
to a family of C's, From the figure, for C = C,, е = Or, at 


МЕ прах“ The entire regenerator ‘temperature distribution is deter- 


mined by following the curve for © = C, from w = O tow š Полу» ШІ 


like manner one can generate curves of m versus w from equation 


дара 





(4,39). In this case the mass flow distribution can be determined by 


following the appropriate C curve from w Z O to v Z « 


па TT addition 


to specifying a range of C's, the regenerator matrix must be known as 
well as the working gas, so that values of n and Y can be assigned in 


equation (A.57). 


FIGURE XIV 


ILLUSTRATION OF COMPUTER SOLUTIONS FOR REGENERATOR 


max 


Spee 





"pulse tube" 


If equation (A.52) is solved for зе 
aw 














Ф Х-1 
СЕТЕ 
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— (4.58) 
2 
m 

Also 

dm _ 1 Фс 

аи а | ае е " Ao 

= ћи qe (4.55) 


Again, the simplest method of solution of these equations is by 
computer, Similar to the case for the regenerator, the simplest approach 
is to treat Ф and C as parameters, Thus, for a given Ф= Q. > 
C is assigned a range of values; then a new value is chosen for 中 and 
C is assigned another range of values. The result is a family of curves 
of C for a given value of ф y as is shown in Figure XV. The curves are 
used in exactly the same manner as that described for the regenerator. 
To use the curves one must be able to determine Ф in advance, From 
equation (4.53) it is apparent that to determine (D , the following 
information must be known: 

(1) the working gas, 

(2) the geometry of the pulse tube, 

(3) the cycling rate (See Chapter IV), 

(4) the temperature of the cooling mediun, 

(5) the average pressure at which the device will be operated 


(defined by equation (A,51)). 


Details of the computer program for the regenerator equations can 


be found in (li). Details of the computer program for the pulse tube 
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cquations can be found in Appendix B, Computer solutions for the гереп- 


erator and pulse tube can be found in Appendix C, 


FIGURE XV 
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APPENDIX B 


General 


Rea (h) wrote a program in Fortran to solve the regenerator equa- 
tions (2.29) and (2.30). A similar program was written to solve the 
pulse tube equations (2.33) and (2.3h). Runge's recurrence formula (6) 
was used in the stepwise integration. 

The program was written for X z 1.67. FUNF(B,R) is equation (2.33), 
where B corresponds to the mass flow, m, and R corresponds to the dimen- 
sionless temperature ratio, ©. The other symbols used in the program 
ere: 

PAR - C 


Ф 
e 


SPAR 


TEMP = T/T, 

Y = m = m/m, 

SLOPE  сФ/4и 

DS an Runge recurrence formula 

E = O” in Runge recurrence formula 

H = integration step size + 0.01 

K = number of integration steps 

KK = number of values of C for each value of 中 

KKK = number of values of Ф 

The first data card contains the initial values of € and m, the ` 


Step size, the number of integration steps, the number of values of C 


for each value of ® , and the number of values of Q . The remaining 
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data cards contain the values of @ and C for which results are desired. 


Comuter Program 
— es A A RO 


Жж жож ж 


172, 


20 
25 


10 
12 
13 
um 


` 
“ 


1.0 


XM 

LABEL 

LISTS 

FORTRAN 
FUNF(B,R)Z(PAR« (1. *SPAR/R )- (0,1:*B) )/(B*B) 

DIMENSION A( JG РР И) SOPEC ЈУ) 
READ 10 TEMP (1), Y(1),H, K, KK, KKK 

DO 25 Lol, KKK 

READ 12,A(L) 
SPARZA(L) 
BOS T= ЕК 
READ 12 с(т) 
PARZ =c(T) 
00.20 J=1;K 

N=¥(5)+(He(1. /TEMP(J ))%*(1.-( (SPAR#PAR)/(TEMP(J )#¥(J))))/2e) 
EZTEMP(J )«( EXFUNF (Y(J ) , TEMP(J ) )/G. 
TEMP(J41)ZTEMP(J ) «CHXFUNF (D, E) ) 
Y(J+1)= Y(J)+(H*(1, /E)#(1,-((SPARzPAR)/(EzD)))) 
W(J+1)=HPLOATE(J) 
SLOPE(J+1 )=FUNF(Y(J+1),TEMP(J+1)) 
PRINT 13,H,PAR,SPAR 
PRINT 15, (W(J41), TEMP(J+1),Y(J+1.),SLOPE(J+1),J=20,K,20) 
CALL EXIT 
FORMAT (3F12.6,1h,12,12) 
FORMAT (1712 .6) 
FORMAT (3F12.6) - 
FORMAT (4F15.5) 

END 

DATA 

170 SO 600 710 


(OTHER DATA FOLLOWS) 
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APPENDIX C 


THEORETICAL AXIAL TEMPERATURE 


Q-PART MODEL. 
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APPENDIX D 


DESCRIPTION OF THE TEST APPARATUS 


The test apparatus consisted of an inlet/exhaust manifold, Pulse 
Tube, heater, heat exchanger, and vacuum chamber. A schematic of the 
apparatus is showm in Figure IV of Chapter IV. 

The inlet/exhaust manifold is shorn in Figure XXIX. The manifold 
was constructed of standard =" copper-tubing, using both soft soldered 
and flared fitting connections. Four Asco solenoid-valves were used to 
control the flow, two on the supply side and two on the exhaust side of 
“һе manifold, These valves had 3/32" orifice. The solenoid valves were 
cycled by a cam operated micro-switch arrangement. The cams were belt- 
driven by a small, electric motor through a four-step pulley which 
ellowed cycling rates of 22.4, Ilı.7, 82.3, and 158.0 reversals per minute. 
"vo 3/8" Hoke needle-valves, one upstream of the inlet solenoid-valves 
and one downstream of the exit solenoid-valves, provided fine control 
of the gas flow rate. The check valve/manual valve loops allowed for 
reading meximum and minimum pressures without any pressure gauge 
fluctuations The maximum pressure was read on a lj"-Marsh, 0-300 psi 
gauge. The minimum pressure was read on 2 3"=#Marsh, O-100 psi gauge. 
The copper~tubing of the manifold mated to a 1 0.D. by 23"~long, hollow, 
brass cylinder. The lower end of this cylinder was machined to allow 
an "O" ring to be inserted, This end of the cylinder was fitted into 
the open-end of the regenerator. Several short lengths of tubing were 
fitted inside the brass cylincer to act as flow smoothers. 


The Pulse Tube drawn to scale is shown in Figure XXX, It consisted 
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FIGURE XXX 


SCALE DRAWING OF PULSE TUBE 
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(10) Copper-Constantan 
Thermocouples 
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SCALE: 4! = 1,0" 








of a regenerator a void tube (pulse tube), and a short length of connect- 
ing tubing. The regenerator was constructed from a i2"-long, siainless- 
steel tube, 1" 0.D. by 0.020"-wall-thiclmess packed with 0,050"~diameter 
lead shot. The pulse tube consisted of a 12"~lonzg, stainless-steel tube, 
5/81 0,2, by 0,016"-wall-thickness. The regenerator and pulse tube were 
connected at the bottom by а 2" 0,0, by 0,020"-wall~thickness, stainless- 
steel tube. On the pulse tube end of the connecting tubing was a 1" 1.0. 
by l"-long, stainless-steel cylinder, which was filled with fine-mesh, 
bronze screens to serve as a flow smoothing device. Several layers of 
this screening were also placed at the bottom of the regenerator for the 
Same purpose. The stainless-steel components were heli-arc welded as 
shown in Figure XXX. The top ends of the regenerator and pulse tube 
were soft soldered 3" epart to a brass, vacuum-jacket, cover-plate, which 
was 9" OD. by 3/8"-thick, Eleven copver-constantan thermocouples were 
installed along the regenerator and pulse tube at locations indicated in 
Figure ХХХ. All the thermocouples were led out of the vacuum chamber 
shroush two l"-diameter Latronic ceramic seals, which were soldered into 
pre-drilled holes in the cover plate, The thermocouple outputs were 
recorded by a Leeds and Northrup multipoint recorder. The pressure- 
time history within the Pulse Tube was measured by diaphragm-type 
pressure transducers mounted at the top of the regenerator and pulse 
tube. An oscilloscove was used to indicate the output of the pressure 
transducers e 

The heater was constructed as follows: the 4! stainless-steel 
tubing connecting the regenerator and pulse tube was coated with a thin 
layer of "Sauereisen Insa-Lute Adhesive Cement", a liquid-porcelain type 


of material. After the cement had hardened, approximately two meters of 


S 





2 high resistance wire were wrapped on top of it. The heater wire was 
connected to two 32-gauge, copper, lead-in wires. Then the entire 
heater assembly was coated again with a second layer of cement. The 
lead-in wires were lead through a Latronic ceramic seal installed in 
the vacuum» jacket cover plate. The total resistance of the heater cir- 
cuit wos 121 ohms. The heater circuit external to the vacuum jacket 
consisted of a 0-1Щ0 volt variac for fine voltage control, A milli- 
ampmeter, voltmeter, and wattmeter were installed in the external cir- 
cuit for power measurements. 

The heat exchanger at the top end of the pulse tube was constructed 
as shown in Figure XXXI. Thin-wall, copper-nickel tubing 1/8" 0,0. ма5 
used to promote good heat transfer. “he working gas was funneled into 
the tubes, which were approximately 34" in length and which were dead- 
ended in a brass tube sheet approximately 1" in diameter by i'"-thick. 
ihe entire tube assembly was surrounded by a water jacket, which was 


attached to the cover plate. Cooling water was pumped through the water 


jacket at the rate of 0.5 gallons per minute. 
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FIGURE XXXI 


HEAT EXCHANGER 
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